ON THE UNIQUENESS OF GIBBS STATES IN THE 
PIROGOV-SINAI THEORY* 



J.L. Lebowitz", and A.E. Mazel ab 

Abstract. We prove that, for low-temperature systems considered in the Pirogov-Sinai 
theory, uniqueness in the class of translation-periodic Gibbs states implies global unique- 
ness, i.e. the absence of any non-periodic Gibbs state. The approach to this infinite volume 
state is exponentially fast. 
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Dedicated to the memory of Roland Dobrushin 
1. Introduction 

The problem of uniqueness of Gibbs states was one of R.L.Dobrushin's favorite sub- 
jects in which he obtained many classical results. In particular when two or more transla- 
tion-periodic states coexist, it is natural to ask whether there might also exist other, non 
translation-periodic, Gibbs states, which approach asymptotically, in different spatial di- 
rections, the translation periodic ones. The affirmative answer to this question was given 
by R.L. Dobrushin with his famous construction of such states for the Ising model, using 
± boundary conditions, in three and higher dimensions [D]. Here we consider the opposite 
situation: we will prove that in the regions of the low-temperature phase diagram where 
there is a unique translation-periodic Gibbs state one actually has global uniqueness of the 
limit Gibbs state. Moreover we show that, uniformly in boundary conditions, the finite 
volume probability of any local event tends to its infinite volume limit value exponentially 
fast in the diameter of the domain. 

The first results concerning this problem in the framework of the Pirogov-Sinai the- 
ory [PS] were obtained by R.L. Dobrushin and E.A.Pecherski in [DP]. The Pirogov-Sinai 
theory describes the low-temperature phase diagram of a wide class of spin lattice models, 
i.e. it determines all their translation-periodic limit Gibbs states [PS], [Z]. The results of 
[DP], corrected and extended in [Sh], imply that, for any values of parameters at which the 
model has a unique ground state, the Gibbs state is unique for sufficiently small tempera- 
tures. But the closer these parameters are to the points with non unique ground state the 
smaller the temperature for which uniqueness of the Gibbs state is given by this method. 
Independently, an alternative method leading to similar results was developed in [Ml, 2]. 

The main difficulty in establishing the results of this type is due to the necessity of 
having sufficiently detailed knowledge of the partition function in a finite domain with an 



a Department of Mathematics and Physics, Rutgers University, New Brunswick, NJ 
08903, USA 

b International Institute of Earthquake Prediction Theory and Mathematical Geo- 
physics, Russian Academy of Sciences, Moscow 113556, RUSSIA 
* Work supported by NSF grant DMR 92-13424 



1 



arbitrary boundary condition. This usually requires a detailed analysis of the geometry 
of the so called boundary /overproduced by such a boundary condition (see [Ml, 2], [Sh]). 
Here we develop a new simplified approach to the problem. The simplification is achieved 
by transforming questions concerning the finite volume Gibbs measure with an arbitrary 
boundary conditions into questions concerning the distribution with a stable (in the sense 
of [Z]) boundary condition. The latter can be easily investigated by means of the polymer 
expansion constructed for it in the Pirogov-Sinai theory. This also allows the extension of 
the uniqueness results from systems with a unique ground state to the case with several 
ground states but unique stable ground state (see [Z]). 

Since the publication of the paper [PS] about twenty years ago the Pirogov-Sinai 
theory was extended in different directions. For a good exposition of the initial theory 
we refer the reader to [Si] and [SI]. Some of the generalizations can be found in [BKL], 
[BS], [DS], [DZ], [HKZ] and [P]. Below we present our results in the standard settings of 
[PS] and [Z]: a, finite spin space with a translation-periodic finite potential of finite range, 
a finite degeneracy of the ground state and a stability of the ground states expressed via 
the so called Peierls or Gertzik- Pirogov-Sinai condition. The extension to other cases is 
straightforward. Our method also works for unbounded spins, see [LM]. 

2. Models and Results. 

The models are defined on some lattice, which for the sake of simplicity we take to be 
the d-dimensional (d > 2) cubic lattice Z d . The spin variable a x associated with the lattice 
site x takes values from the finite set S = {1, 2, . . . , \S\}. The energy of the configuration 
a E S z is given by the formal Hamiltonian 

H (a)= J2 U ^a). (1) 

iCZ d , diam A<r 

Here a a £ S A is a configuration in A C Z d and the potential Uq{<ja) '■ S A i— > R, satisfies 
Uo(o~a) = Uo(o~A+y) for any y belonging to some subgroup of Z d of finite index and the sum 
is extended over subsets A of Z d with a diameter not exceeding r. Accordingly for a finite 
domain 7cZ d , with the boundary condition a yc given on its complement V c = Z d \ V, 
the conditional Hamiltonian is 

H (a v \a vc )= Yl Moa), (2) 

Ar\V^0, diam A<r 

where a a = ctadv + o~adv c for A(~)V C ^ 0, i.e. the spin at site x is equal to a x for x E An V 
and a x for x E A fl V c . 

A ground state of (1) is a configuration a in Z d whose energy cannot be lowered by 
changing a in some local region. We assume that (1) has a finite number of translation- 
periodic (i.e. invariant under the action of some subgroup of Z d of finite index) ground 
states. By a standard trick of partitioning the lattice into disjoint cubes Q(y) centered 
at y E qZ d with an appropriate q and enlarging the spin space from S to one can 
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transform the model above into a model on qZ d with a translation-invariant potential and 
only translation-invariant or non-periodic ground states. Hence, without loss of generality, 
we assume translation-invariance instead of translation-periodicity and we permute the 
spin so that the ground states of the model will be . . . , a*™) with = k for any 
x G Z d . Taking q > r one obtains a model with nearest neighbor and next nearest neighbor 
(diagonal) interaction, i.e. the potential is not vanishing only on lattice cubes Qi of linear 
size 1, containing 2 d sites. 

Given a configuration a in Z d we say that site x is in the k-th phase if this configuration 
coincides with inside the lattice cube Q2(x) of linear size 2 centered at x. Every 
connected component of sites not in one of the phases is called a contour of the configuration 

a. It is clear that for a = a xr + a^ k ) contours are connected subsets of V which we denote 

v v 

by 7i(cr), . . . , 7;(<r). The important observation is that the excess energy of a configuration 
a with respect to the energy of the ground state is concentrated along the contours 
of a. More precisely, 

Ho (a v | o$ ) - H (*W \ a ) = £ H {% (a) ) , (3) 

i 

where 

H Q {li{a))= {U (<T Q1 ) - U (o- { £)) (4) 

Ql_: QiCji(a) 

and the sum is taken over the unit lattice cubes Qi, containing 2 d sites. The Peierls 
condition is 

#o(7^))>r|7^)|, (5) 
where r > is an absolute constant and |7«(cr)| denotes the number of sites in 7i(cr). 

Consider now a family of Hamiltonians H n {o) = Uni&A), n= 1 5 • • • , ™~ 

AcZ d , diam A<r 

1 satisfying the same conditions as Hq with the same or smaller set of translation-periodic 
ground states. For A = (Ai, . . . , A m _i) belonging to a neighborhood of the origin in R m_1 
define a perturbed formal Hamiltonian 

m — l 

H = H +J2 ^H n . (6) 

n=l 

Here X n H n play the role of generalized magnetic fields removing the degeneracy of the 
ground state. The finite volume Gibbs distribution is 



exp [-(3H(a v \a vc )] 
*{V\a vc ) 



^V^V^ a v) = ^nr^ ^ ' ( 7 ) 



where (3 > is the inverse temperature and liy^yAoy) is the probability of the event 
that the configuration in V is o y , given o yc . Here the conditional Hamiltonian is 



m— 1 



H(a v \a vc )= ^( a Qi)i U(-)=J2 U n(-) (8) 

Q 1 :Q 1 r\V^<D n=0 
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and the partition function is 



Z(V\a yc ) = ]Texp [-(3H(a v \a vc )} (9) 
V 



The notion of a stable ground state was introduced in [Z] (see also the next section) 
and it is crucial for the Pirogov-Sinai theory because of the following theorem. 

Theorem [PS] , [Z] . Consider a Hamiltonian H of the form (6) satisfying all the conditions 
above. Then for f3 large enough, (3 > (3o(X), every stable ground state a^ k ' generates a 
translation-invariant Gibbs state 

/*<*>(•)= lim /* («(•)• (10) 

These Gibbs states are different for different k and they are the only translation-periodic 
Gibbs states of the system. 

Remark. Given H n , n = 0, . . . , m — 1 there exists sufficiently small A such that for 
max |A n | < Ao the quantity (3$ becomes independent on A. 

l<n<m— 1 

An obvious corollary of the above theorem is 

Corollary. If there is only a single stable ground state, say a^ 1 ', then for (3 > (3o(\) there 
is a unique translation-periodic Gibbs state 

^(■)= lim u (11) 

V^Z d V'"y c 



Our extension of this result is given by 

Theorem. Consider a Hamiltonian H of the form (6) satisfying all the conditions above. 
Suppose that given A and j3 > Po(X) there exist a single stable ground state, say a^K 
Then the Gibbs state /z(-) = A* (') JS unique. Moreover, for finite A C V, such that 
dist(A, V c ) > 2d diamA + Ci(r, (3, A, d), any configuration a a and any boundary condition 
o~y C , one has 

where C U C 2 > 0. 

Remark. In contrast with [DP], [Sh] and [Ml, 2] the theorem above treats the situa- 
tion when there are several ground states with only one of them being stable. Moreover, 
the result is true for all sufficiently low temperatures not depending on how close the pa- 
rameters are to the points with non unique ground state. If some of the conditions of the 
Theorem are violated the statement can be wrong. 

The simplest counterexample can be constructed from the Ising model in d = 3. 
It is well-known that at low temperatures this model contains precisely two translation- 
invariant Gibbs states taken into each other by ± symmetry and infinitely many non 



< exp [-C 2 (r, (3, A, d) dist (A, V c )] , (12) 
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translation-invariant Gibbs states, i.e. the Dobrushin states mentioned earlier. Identifying 
configurations taken into each other by ± symmetry one obtains a model with unique 
translation-invariant Gibbs state and infinitely many non translation-invariant ones. The 
condition of the Theorem which is not true for this factorized model is the finiteness of 
the potential: the model contains a hard-core constraint. 

Formally speaking, one can consider a model with spin variables still taking values ±1 
but assigned to bonds of the lattice. Then the Hamiltonian is the sum of the spins over all 
lattice bonds multiplied by a negative coupling constant. The hard-core constraint says 
that the product of spins along any lattice plaquette is 1. 

Another example based on a gauge model can be found in [B]. For this model the 
spin space is finite, the potential is finite and of finite radius but the Peierls condition is 
violated and the number of ground states is infinite. On the other hand, after a proper 
factorization, this system can be transformed into a model with a hard-core restriction 
similar to that discussed above. 

It is known [LML] that for systems satisfying FKG inequalities the uniqueness of the 
translation-invariant Gibbs state implies uniqueness of all Gibbs states at any temperature. 
This makes it tempting to assume that the conclusions of our theorem hold beyond the 
low temperature region covered by the Pirogov-Sinai theory. We are not aware of any 
counterexample. 

3. Proof of the Theorem 

In this section we assume that the reader is familiar with the Pirogov-Sinai theory and 
we only list the appropriate notation. Then we quote some necessary results and proceed 
to the proof of the theorem. 

Preliminaries. A contour is a pair 7 = (7, a^) consisting of the support 7 and the 
configuration in it. The components of the interior of the contour 7 are denoted by 
Intj^f and the exterior of 7 is denoted by Ext'y. The family {7, Intj^, Ext^f} is a partition 
of Z d . The configuration can be uniquely extended to the configuration a' in Z d taking 
constant values -^(7) and £'(7) on the connected components of 7°. Generally these 
values are different for different components. The contour 7 is said to be from the phase k 
if £'(7) = k. The energy of the contour 7 is 





(13) 



Qi: QiC7(cr) 



The statistical weight of 7 is 



w(7) = ex P(-/^(7)) 



(14) 



and satisfies 



< tuft) < e-^L 



(15) 



5 



The renormalized statistical weight of the contour is 



„„ , , , rr 5(/r.t; 7 |/j(7)) 
W(7) = t " (7> n 5(/raf ; 7 | E(7) ) 



where for any A C Z d we denote 

A* = {x G A\ x is not adjacent to A c }. 
The contour 7 is stable if 

1 



(16) 



(17) 



(18) 



and the ground state a^ k ' is stable if all contours 7 with E(p() = k are stable. It is known 
(see [Z]) that at least one of the ground states is stable. Because of (18) for any x G Z d , 
N > 1 and (3 large enough 

E Wirt) < e~ c > N , (19) 

7: (Ext-y) c 3x, |7|>7V 

where C3 = C3(r, /3, d) is positive and monotone increasing in r and (3. In particular 



E ^(7) < C4, 

7: {Ext~i) c 3x 



(20) 



where C4 = e C3 . For the stable ground state cr( fe ) the corresponding partition function 
can be represented as 



~(V\k) = e-^v^ 



e n^(^). 

[7i] s eV, B([ 7i ]»)=fe i 



(21) 



where the sum is taken over all collections of contours [7$] such that 7$ are disjoint, i?(7i) = 
for all % and 7i C y for all i. 

Representation (21) and estimate (18) allow to write an absolutely convergent polymer 
expansion 

\og~(V\k) = -(1H(v ( y ) W { y})+ E W(* (k) ), (22) 

where the sum is taken over so called polymers of the phase k belonging to the domain 
V. By definition a polymer = (7$) is a collection of, not necessarily different, contours 
7i of the phase such that VJi^i is connected. The statistical weight W(tt^) is uniquely 
defined via W{^i) and satisfies the estimate (see [Se]) 



\W(-n {k) ) \ < exp 



\pr - 6d) E 



Til 



(23) 
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implying 

\W(tt^)\ < e~ C:iN . (24) 

Denote by ({7i}> ext) the probability of the event that all contours of the collection 
{7^} are external ones inside V. By the construction 

tfHHhexQKjlWbi). (25) 

i 

From the polymer expansion (22) and estimate (24) it is not hard to conclude that for {7^} 
with dist(Ui7 l7 V c ) > I U,7i| 

^yHHhext) -// fc >({ 7i },ex*)| < / /( fc )({ 7i },ea;t)| U i7i |exp[-C 5 dist (U l7l , V% (26) 

where C5 = C^(T,P,d) is positive and monotone increasing in r and (3. For any A G V, 
dist(A, F c ) > diamA, and any cr^ estimate (26) implies in a standard way that 

ltf> (a A ) - {a A )\< exp [-C 6 dist (A, F c )] , (27) 

where again C§ = Cq(t, [3, d) is positive and monotone increasing in r and /?. 

For an arbitrary boundary condition a yc the probability \i y ^ yc {? a) depends only on 
a , where 

ay = G F c : x is adjacent to V}, (28) 

and we freely use the notation fi v ^ (&a)- From now on we suppose that is the 
only stable ground state of H and denote by fi y (-) the Gibbs distributions with the stable 
boundary condition cr yc - 

For a domain V with the boundary condition a^] fix / < L sites belonging to dV 
and consider a collection { 7 i} e G V of external contours touching <9F at one of these sites. 
A smaller domain V = Uj L)j Intj^i has a natural boundary condition induced by 
UjU^. Given M > J7 | 7 j| denote by £{ 7 .}e ^ the event that the total number of sites in 
all adjacent to dV connected components of sites not in the 1-st phase is not less than 
M — J2 i | 7 i|. According to the Theorem of Section 3.2 in [Z] 

Vv'Kv fa?*}'.") ^ e- C7 ^' A ' d) ( M -Sj^l )+C8 ^' A ' d ^jT,| (29) 

The positive constants C7 and Cs tend to as (3 — > 00 or (/?, A) approaches the manifold 
on which a*- 1 -* is not the only stable ground state. For different { 7 i} e the events £{ 7 .}e M 
are disjoint and for their union £v,m,l = ^{-yi} e ev^{-yi} e ,m one has the estimate 

V>y (£V,M,L) = Yl Vy',*' ev , {^{^YM) Vy {{liY) 

{nYev 
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< J- e -C 7 (M-E,|7,|)+C 8 |7,|JJ^ (7j) 
< e -C 7 M -Q e _(i T/3 _c 7 -C 8 )| 7l | 

< e - c ' M (l + C 4 ) L 

Finally observe that for any A C.V and any a a 

V v (a A )e- c » L ^ < fJL v , &vc ((T A ) < V v {° A )e C » L ^ v \ (31) 
where L(asv) is the number of sites x G dV with oqv 7^ 1 an d 

C 9 = 2 d /3max|t/( ( 7Q 1 )|. (32) 



Proof. We are now ready to prove the theorem. Take an integer N > and suppose that 

V contains a cube Qgn with sides of length 6N centered at the origin. From now on all 
cubes are assumed to be centered at the origin. Let Qn', N' > 6N be the maximal cube 
contained in V. Denote d'V = dQ^i fl dV. First we consider boundary conditions a gv 

which coincide with on dV \d'V and differ from on d'V by at most y/~N lattice 
sites. 

Given a y denote by 0(cr ) the union of the connected components of the set {x G 

V : x is not in the 1-st phase } adjacent to {x G dV : a x ^ 1}. This set is called 
the boundary layer of o~ y and we denote by Qi(o~ y ) its connected components. Introduce 
the event £ = {a y : 0(cr y ) fl Q 4 n ^ 0}- By construction for a y G So every fl^tr ) 
touches dV at some site x G dV with a x ^ 1 and there exists at least one component 
intersecting Q 4 n- Without loss of generality we suppose that it is Qi(a y ). This leads to 
the estimate 



^ yc (£o)<e c ^ N ^ y (£o) 

<e c ^VNe- CsN (l + C 4 )^ 

< e ~ CloN . (33) 

In the first inequality of (33) we used (31) reducing the problem to the calculation 
for the stable boundary condition The second inequality comes in a standard way 

from the cluster expansion for fi y (-). Indeed, in the domain V with the stable boundary 
condition a^ y } every component Oj contains an external contour 7^ such that E{p/i) = 1, 
7i C fij and C (Ext(^i)) c . One may simply say that 7$ is the external boundary of 
Cli and clearly 7, touches dV. If fii intersects Q 4 n then 71 intersects or encloses Qan- 
The number of possibilities to chose the site x G dV, a x 7^ 1 at which 71 touches dV 
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does not exceed V N which produces the factor VAT in the estimate. The next factor 
estimates the sum of the statistical weights of all possible 71 touching this site. It is based 
on (19) and takes into account the fact that the diameter of 71, and hence I71I, is not 
less than N. The constant C4 (see (20)) estimates the sum of the statistical weights of 

all possible 7^ touching given lattice site and (1 + C^f^ estimates the statistical weight 
of all possibilities to choose {7^ i 7^ 1}. The whole estimate uses (25) and the fact 
that fi v ({7^}, ext) is the upper bound for the sum of //^-probabilities of boundary layers 
O = {Oi} having 7^ as the boundary of O^. The third inequality in (33) is trivial for 
C10 = 0.5C 3 and N > 4(C 9 + 1 + C 4 ) 2 C 3 " 2 . 

Denote by Sq the complement of Sq. If o y G Eq then (V \ Q(a y ))* D Qan-2 (see (17) 
for the definition of (■)*)• It is not hard to see that the configuration a y G S% equals 1 on 
the boundary of (V \ Q(a v ))*. Now fix A C Q2N and a a- In view of (27) one has 



< e 



-C 6 N 



(34) 



This gives us 



< 



+ 



V' a V 



(a A \£ ) - fx(a A ) 



V v ,a, /c {<TA\£o) -K<7A) 



< e 

< e 



"V> a V 
-C W N 

-CnJV 



fj, 



y>°V' 



:(S>) 

(£0) 



+ e 



-C 6 N 



(35) 

where Cn = 0.5min(Ci , C 6 ) and N > log2/Cu. 

To extend (35) to the wider class of boundary conditions we suppose that V D Qsn 
and Qn'i N' > 8N is the maximal cube contained in V. Now we consider boundary 
conditions a which coincide with on dV \ d'V and differ from crW on d'V by 

dV x J 

at most (VN) 2 lattice sites. Denote Aj = QsN-2i+2 \ QsN-2i and let Q^\a v ) be a 
union of the connected components of the set {x G V \ Qsn-2i '■ &x 7^ 1} adjacent to 
{x G dV : a x 7^ 1}, i = 1, . . . , N. Introduce disjoint events 

S x = {a v : |n^)(0 n A x | <v / iV}, 



Si = {a v : \n w (a v ) n Ai| > Vn, . . . , |ft (i - 1} (<v) n a<_i| > Vn, 



v 



V 



n {l) (a v )n A, 1 < Vn} 



and 



(36) 



If cr y G £ c than the boundary contour n(cr y ) contains at least NyN sites. Hence (30) 
implies the following estimate for the probability of £ c 

*V,* V c(*c) < e C9(vW)2 /i v (£ c ) 

< e C 9 (v / iV) 2 e -C , 7 Afv / lV+C4(v / iV) 2 

< e" Cl2iVvW , (37) 
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where C 12 = 0.5C 7 and N > 4(C 9 + C 4 ) 2 Cf 2 . 

For a £ Si consider the volume Vi = (V \ ftW (a ))* U Qsn-2% with the boundary 
condition o yc + cf v \ v .. By construction the number of sites x G dVi with a x 7^ 1 is less 



than v N and one can apply (35) to obtain the bound 



f* v ,<t vc (<TA\£i) - fi(<TA) 

Joining (37) and (38) we conclude 



< e 



-CuN 



(38) 



»v> a v 



< e -CuN + e -C 12 Ny/N 



< 2e 



-CuAT 



(39) 

where AT > (C n /Ci 2 ) 2 . Expression (39) is a version of (35) which is weaker by the factor 
2 in the RHS but is applicable to the wider class of boundary conditions containing (v^V) 2 
unstable sites instead of VW for (35). 

The argument leading from (35) to (39) can be iterated several times. The first 
iteration treats the following situation. Suppose that V D Qiotv and Qn>, N' > ION is 
the maximal cube contained in V. Consider boundary conditions a which coincide with 

cr^ 1 ) on dV \ d'V and differ from on d'V by at most (VW) 3 lattice sites. Then the 
analogue of (39) is 



< 3e~ 



(40) 



Similarly after 2d iterations one obtains that for any V D Q 2c in with the boundary condi- 
tion a containing not more than {\fN) 2d unstable sites 



fj, 



v>°~v c 



{a a) - v((7a) 



< 2de~ c ^ N . 



(41) 



Now set Ci = 2d max (4(C 9 + 1 + C 4 ) 2 C^ 2 , log2/Cu, 4(C 9 + C 4 ) 2 Cf 2 , (Cu/C 12 ) 2 ) and 
for any A and a a consider a cube Q 2 l with L > C\ and dist(A, Q 2L ) > (1 — l/d)L. Taking 
N = L/d and C 2 = C n /2d for V = Q 2L one obtains (12) from (41). For any V D Q 2L 
with dV fl 8Q 2 l 7^ we have 



< exp [— C 2 dist (A, F c )] , (42) 
which finishes the proof of the Theorem. 

Acknowledgments: We are greatly indebted to F.Cesi, F.Martinelli and S.Shlosman 
for many very helpful discussions. 
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